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DELTA7 :  A  COMPUTER  MODEL  OF  DELTA  GROWTH 


Introduction  and  Summary 


\  V-.  --'" 


De1ta7  Is  an  Interactive  Fortran  IV  Computer  Program  written 

\ 

for  the  PDP-11  mini -computer.  It  accepts  initial  beach  and  river 
conditions,  as  well  as  a  varying  set  of  deepwater  wave  conditions, 
and  produces  a  data  set  that  outlines  a  planview  of  the  coastline 
at  specified  time  intervals.  The  coastline  can  subsequently  be 
graphed  using  the  Calcomp  1037  Drum  Plotter  with  the  XPLOT  instruc¬ 
tion  of  MAGIC. 


\ 


Two  main  background  sources  in  the  literature  were  found  valuable 
in  developing  Delta7.  Komar  (1973)  studied  the  problem  of  the  growth 
of  a  river  delta,  and  Price,  Tomlinson  and  Wilson  (1972)  used  a 
computer  to  simulate  coastal  processes  around  a  jetty.  Kim-E  (1978) 
attempted  to  duplicate  Komar's  data,  as  well  as  to  add  an  approximation 
for  refraction.  Those  programs  and  the  discussion  presented  were  used 
as  a  starting  point  in  developing  Delta7.  Delta7  is  then  a  continuation 
of  Kim-E's  extension  of  Komar's  original  work. 

The  main  objective  in  developing  Delta7  was  to  determine  if  varying 
the  order  of  attacking  waves  would  significantly  affect  the  coastline 
development.  It  was  strongly  suspected  that  it  would,  and  therefore 
that  the  practicality  of  a  computer  simulation  would  be  limited,  due 
to  the  necessity  of  generating  ordered  wave  data  for  the  future!  An 
implied  objective  was  to  verify  Delta7  by  duplicating  Komar's  data  on 
one  year  delta  growth  for  a  constant  wave  of  P  =  5.0  x  10s  ergs/cm-sec, 
at  a  10  degree  angle  (see  Komar  (1973)).  Also  implied  was  the  requirement 
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to  find  a  realistic  wave  data  set  for  input.  To  this  end,  564  days 
worth  of  the  most  frequently  occurring  waves  from  Wave  Station  5  off 
the  Southern  California  Coast  (see  Meteorology  International  (1977)) 
were  extracted  and  put  into  a  file  called  WAVE1.DAT,  which  is  listed 
In  Its  entirety  in  Appendix  F.  This  is  somewhat  of  an  expedient  data 
base,  since  It  excludes  all  the  extreme  waves  outside  of  the  40  degree 
fan  within  which  aQ  most  often  occurs,  and  within  that  fan  waves  with  a 
height  greater  than  2.5  meters.  It  was  felt,  however,  that  the  main 
objective  could  be  reached  without  adding  the  entire  spectrum  into  the 
input,  and  the  results  bear  this  out.  After  the  data  set  WAVE1.DAT  was 
created  it  was  "fl ipped"— that  is,  its  order  was  reversed.  (The  original 
data  being  in  order  of  increasing  wave  size.)  It  was  then  shuffled  three 
times  to  produce  three  data  sets  in  random  order.  The  different  versions 
of  WAVE1.DAT  were  then  used  as  input  to  Delta7  producing  several  different 
coastlines,  all  of  which  had  been  generated  by  the  same  waves,  arriving 
in  different  order.  Refer  to  Appendix  F  for  Figures  15,  16,  19  which 
show  the  ordered  data  graphically  and  a  representation  of  the  different 
random  curves  over  a  selected  period  of  time. 

Analysis 

The  simulation  uses  a  model  which  divides  the  beach/coast  into  an 
infinite  series  of  cells,  which  are  parallelepipeds  of  equal  width  Ax, 
equal  depth  d,  which  may  vary  linearly  if  the  beach  has  a  slope,  and  a 
variable  height  y,  which  Is  where  accretion  or  erosion  is  reflected.  It 
Is  assumed  at  time  t  *  0,  that  all  y  =  0,  that  is  we  are  starting  with 
a  straight  coastline.  Sand  is  then  introduced  by  the  river  at  a  rate 
of  20,000  cubic  meters  per  day.  Longshore  transport  distributes  the 


sand  along  the  coastline,  building  the  delta,  or  in  some  cases,  washing 
It  away.  It  Is  assumed  that  no  sand  is  lost  between  the  cells,  or 
gained,  other  than  by  longshore  transport  or  from  the  river  sediment. 

We,  therefore,  neglect  wind  transfer,  onshore-offshore  losses,  erosion  of 
larger  particles,  etc.  One  may  then  write  a  conservation  of  sand 
equation  for  each  cell  (see  Komar  (1976)). 


(1) 


where  A V  is  the  change  in  volume  in  cubic  meters.  At  is  the  change  in  time, 
in  days,  and  AS  is  the  rate  of  change  In  sand  volume  for  the  cell  In  cubic 

meters  per  day.  Substituting  an  expression  for  the  volume  of  the  cell,  we 
develop  an  expression  for  the  change  in  the  height  of  the  cell  or  the 
change  in  the  coastline. 


Ay  = 


(S..!  -  S^At 
dAx 


(2) 


where  Ay  represents  the  distance  the  coastline  erodes  or  accretes.  In 
Equation  (2),  we  set  At,  Ax  and  d  constant;  and  therefore  we  need  only 
to  define  S^  _-|  and  S..  to  determine  Ay,  and  thence,  by  iteration,  the 
entire  coastline.  An  empirical  expression,  using  the  longshore  current, 
is  given  by  Komar  (1976).  Starting  with  an  equation  energy  flux  in  the 
alongshore  direction 

1^  =  K(ECn)cos0s1n6  (3) 

where  K  is  an  empirically  determined  constant,  given  by  Komar  (1976)  as 
.77,  ECn  is  an  expression  for  wave  power  in  ergs/cm-sec,  and  e  is  the 
angle  the  breaking  wave  makes  with  the  shoreline.  Further,  we  can  write 
an  expression  for  the  immersed  weight  transport  rate  as: 
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ll  =  (ps  "  p)ga'S  (4) 

where  p$  for  a  typical  beach  sand  was  taken  as  2.65  grams/cm3,  p  was 
taken  as  1.00  grams/cm3  for  water,  g  was  the  constant  at  sea  level 
of  980  cm/sec2,  and  a'  is  another  empirically  determined  constant  for 
the  porosity  of  the  sand,  which  was  set  at  .61.  In  order  to  get  S  in 
units  of  m3/day,  we  must  also  add  a  conversion  for  cm3/sec.  The 
result,  combining  Equations  (3)  and  (4)  is 

S  _  K(ECn)cos9s1n6  ^ 

(p$-p)ga' 

.77(ECn)  cos9sin6  (10~6  m3/cm3) (36,400  sec/day 
.65  gm/cm3)(980  cm/ sec2)  .61 


S  =  6.85  x  1O"5(ECn)cos0sin0  (m3/day)  (7) 

At  this  point,  it  will  be  worthwhile  to  discuss  the  angle  9  in  some 
detail,  since  it  was  the  source  of  some  confusion  during  the  develop¬ 
ment  of  the  model.  Referring  to  Figure  1,  0  is  the  angle  the  breaking 
wave  makes  with  the  coastline.  The  other  angles  are  all  relative  to 
the  x-axis,  which  was  the  coastline  at  t  =  0;  our  starting  assumption. 
aQ  is  the  angle  between  the  crest  of  the  refracted  deepwater  wave  and 
the  x-axis,  and  6  the  angle  the  contour  of  the  coastline  makes  with 
the  x-axis.  Using  small  amplitude  wave  theory,  the  angle  ab  is 
obtained  as 

si  net-  sina. 

_ O  -  _ b  fox 


wherein  the  ratio 
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=  tanh  (kh).  (9) 

uo  D 

combining  Equations  (8)  and  (9)  the  following  equation  results  for 
small  values  of  kh,  i.e.,  for  breaking  conditions: 

sina^  =  (kh)b  sinaQ  (10) 

since  for  small  values  of  x,  tanh  x  =  x.  Finally,  a  relationship  for 
( kh ) b  was  developed  using  the  fact  that  the  power  between  orthogonal s 
remains  constant  when  the  wave  is  refracted: 

2/5  1/5 

(kh)b  -  1.901 (H0/L0)  (cosaQ)  (11) 

The  derivation  for  Equation  (11)  is  presented  in  Appendix  E.  In 

Equation  (11),  HQ,  LQ  and  aQ  are  respectively  the  deepwater  wave 

height,  wavelength  and  angles  with  the  x-axis.  Using  Equations  (10) 

and  (11),  we  may  then  express  the  angle  of  the  breaking  wave  entirely 

In  terms  of  Its  deepwater  parameters  HQ,  LQ  and  aQ: 

r  2/5  1/5 

ab  =  arc  [sin  1.901  (HQ/L0)  (cosaQ)  sinaQJ  . 

Combining  ab  and  8  we  get  6,  which  Is  the  key  element  in  Equation  (7). 
The  other  part  of  Equation  (7)  can  also  be  written  in  terms  of  the 
deepwater  parameters  In  the  following  expression  from  Kim-E  (1978): 

P  =  ECn  =  9.55  x  107HQ2T  (12) 

where  T  is  the  wave  period  in  seconds  and  Hq  as  before  is  the  deep  water 
significant  wave  height,  in  meters. 


Presentation  and  Discussion  of  Results 


The  actual  program  was  evolved  through  six  revisions,  Delta7 
being  the  last.  Delta4  was  a  program  that  simulated  one  type  of 
wave  for  a  specific  length  of  time,  and  Delta6  was  similar  to  Delta7 
except  that  it  did  not  refract  the  deepwater  wave.  Although  the 
program  of  Kim-E  (1978)  furnished  a  start  point,  it  was  clear  that 
it  was  Inefficient  because  it  used  an  excessive  amount  of  large  arrays 
and  did  twice  the  amount  of  required  calculations.  Using  instead, 
Komar's  expression  for  AS  =  S.._i  -  S..  (see  Komar  (1976)),  it  was 
possible  to  compute  the  approximation  for  A V  only  once  per  iteration 
instead  of  twice  as  had  been  done  previously.  This  had  the  added 
benefit  of  Insuring  that  sand  was  conserved  between  the  cells,  since 
the  sand  out  of  cell  j  automatically  became  the  sand  into  cell  j+1. 

In  addition,  Kim-E  (1978)  had  used  the  difference  between  two 
cells  to  determine  the  angle  B  for  the  coastline.  The  Delta  series 
uses  a  central  difference  as  suggested  by  Weggel  (paper  undated). 

A  central  difference  uses  the  difference  between  the  j-1  and  j+1 
cells  over  2Ax  to  determine  8  at  the  jth  cell.  In  verifying  an 
example  from  Komar's  (1973)  data,  the  program  developed  in  this 
study  and  denoted  as  Delta4  was  closer  than  Kim-E  as  shown  in  Fig.  5, 

primarily  because  of  the  change  to  a  central  difference,  which  has  a 

moderating  effect  on  the  values  of  B.  Finally,  although  Komar  (1973) 

and  Kim-E  (1978)  use  the  trigonometric  Identity  for  tan(a±b),  in  this 

analysis  the  angles  are  added  after  taking  the  arc  tan;  this  results 
in  simpler  processing,  and  avoids  confusion  in  sign  convention. 

Delta7  Is  broken  down  into  a  main  program  and  two  subroutines 
INPUT  and  SWELL  as  illustrated  in  the  flow  diagrams.  Figs.  2  through  4. 
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Figure  2  Main  program  -  Delta7. 
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Figure  3  Flow  diagram  subroutine  swell  -  Delta) 


Figure  5  Comparison  of  Komar,  Klm-E  and  Delta4  for  verification 


12 


A  list  of  variable  names  is  attached  as  Appendix  A,  and  a  complete 
listing  as  Appendices  B  through  D.  INPUT  accepts  beach  data,  and 
program  utility  values  for  accuracy,  intervals  and  total  length  of 
simulation.  SWELL  reads  deepwater  data  from  file  WAVE1.DAT,  computes 
wave  power  In  10®  ergs/cm-sec,  refracts  the  angle  aQ  to  get  afa,  and 
finally  increments  total  duration  by  the  duration  of  the  new  wave. 

The  main  program  functions  within  two  loops;  one  travels  along 
the  beach  and  iterates  for  every  cell,  and  the  other  iterates  the 
entire  beach  for  each  1/10  day.  Both  the  Increments  of  time  and 
beach  length  could  be  changed  easily.  However,  Komar  (1973)  uses 
1/10  day  and  100  meters  and  it  does  not  appear  as  if  the  Increase  in 
accuracy  could  be  justified  by  either  the  accuracy  of  the  input 
information  or  the  increased  cost,  t.e.,  computation  time,  required. 

The  program  moves  first  updrift  and  then  downdrift  from  the  river 
computing  change  for  each  cell  until  the  cell  is  reached  where  growth 
is  less  than  an  Inputted  limit,  generally  .01.  Time  is  then  incremented 
and  the  process  repeated. 

A  new  wave  is  inputted  and  coastlines  are  written  as  determined 
by  the  values  of  variables  TDUR  and  TOTAL,  respectively.  When  time 
exceeds  the  total  duration,  TDUR,  at  the  end  of  a  wave,  SWELL  is  called. 
The  subroutine  automatically  reads  new  wave  parameters,  converts  them 
to  usable  form  and  increments  TDUR.  Similarly,  when  time  exceeds 
TOTAL,  a  coastline  is  printed  for  display  and  written  on  file  DELDAT. 
Then,  before  going  to  the  next  iteration  TOTAL  is  incremented  by  TINTR, 
the  time  interval  between  desired  prints,  in  our  data  generally  90  days. 

Results  are  in  the  attached  graphs.  The  attempt  to  verify  the 
results  of  Komar  (1973)  is  shown  in  Fig.  5  (graph  1).  There  is  close 
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agreement  with  Komar  (1973),  and  at  least  a  relative  improvement  over 
the  verification  presented  by  Kim-E  (1978).  Wave  Station  5  data  were 
used  unrefracted  in  Graphs  2  and  3,  Figs.  6  and  7.  When  the  waves  are 
ordered  by  significant  height  and  period,  thence  by  power,  a  "radical" 
delta  evolves.  With  decreasing  power  an  Initially  flat  delta  becomes 
sharp  and  pointed  (Graph  2,  Fig.  6),  while  with  increasing  power  a 
pointed  delta  is  "pounded"  flat  (Graph  7,  Fig.  9).  However,  when  the 
same  waves  are  randomly  sorted,  they  result  in  what  can  be  termed  a 
"normal"  delta,  which  varies  only  slightly  from  sort  to  sort.  Since 
nature  randomly  sorts  waves  to  a  certain  extent,  it  may  be  possible  to 
use  this  model  to  predict  long  term  coastline  evolution  to  a  degree 
that  would  be  useful  for  engineering  purposes. 

The  addition  of  the  refraction  equation  to  SWELL  had  only  a  minor 
effect  on  Delta  Growth  as  shown  in  Graphs  6  through  10,  probably  because 
the  range  of  values  for  aQ,  the  deepwater  wave  angle  were  scaled  down  to 
the  range  0  to  10  degrees  In  the  input  set.  However,  as  Komar  (1973) 
Indicates, the  influence  of  wave  angle  is  not  as  pronounced  as  we  might 
have  expected,  at  least  over  the  range  indicated.  It  appears  the 
effect  of  wave  power  is  significantly  more  important.  Figure  13  (Graph  11) 
shows  a  Delta7  run  with  a  WAVE1.DAT  data  set  modified  so  that  the  deep¬ 
water  angle  again  ranges  from  0  to  40.  The  waves  are  refracted,  and 
comparison  with  previous  random  graphs  show  no  major  changes.  Figure 
14  shows  a  series  of  graphs  which  indicated  the  marked  effect  ordered 
wave  power  has  on  growth  at  selected  points  on  the  coast. 


Strict  ordering  of  waves  does  significantly  change  coastline 
evolution  with  this  model.  However,  randomly  sorted  waves  tend  to 
produce  approximately  the  same  delta,  which  indicates  it  may  be  possible 


Delta6  with  a  564  day  simulation  using  waves  In  order  of  decreasing  power.  Unrefracted 
Contours  are  at  90  day  Intervals.  Solid  line  Is  final  shoreline. 


Figure  7  A  564  day  simulation  by  Delta6  using  a  random  file  with  deepwater  angles  between 
0  and  10  degrees.  Graph  2  (Figure  6)  was  the  same  data  but  ordered. 
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Figure  8  A  comparison  of  the  differences  in  Delta7  (refracted)  runs.  The  symbols  represent  the 
outline  of  an  ordered  wave  attack,  with  waves  in  order  of  Increasing  power  as  In  Graph  7 
(Figure  9).  The  other  curves  are  the  same  waves  but  in  three  different  random  sorts. 


Delta7  (refracted)  with  a  564  day  simulation  shown  at  90  day  Intervals  with  an  Increasing 
ordered  wave  attack.  Deepwater  angles  from  0  to  10  degrees. 
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Figure  10  De1ta7  (refracted)  with  a  564  day  simulation  with  a  random  wave  attack.  Deepwater 
angles  from  0  to  10  degrees. 
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Figure  13  Delta7  (refracted)  same  as  Figures  10,  11,  12  but  the  waves'  angles  have  been  Increased 
by  a  factor  of  4.  Deepwater  angles  vary  then  from  0  to  40  degrees. 


Figure  14a  Growth  at  fixed  points  on  the  delta  plotted  against  time.  The  solid  line  is  wave  data  ordered 

to  be  increasing  with  respect  to  power.  The  dashed  lines  are  the  three  different  random  sortings 
of  the  same  data.  The  waves  were  refracted.  Data  is  shown  here  for  100,  500,  and  1000  meters 
from  the  river,  while  Figure  14b  has  points  at  1500,  2000  and  2500  meters. 


Figure  14b  The  waves  were  refracted.  Data  is  shown  here  for  1500,  2000  and  2500  meters  from  the  river,  while 
Figure  14a  has  points  at  100,  500  and  1000  meters  from  the  river. 
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to  realistically  simulate  coastline  evolution.  The  problem  remains, 
however,  to  have  a  valid  input  data  set,  and  this  may  have  to  be  somewhat 
ordered  to  reflect  at  least  winter  and  summer  wave  conditions.  Wave 
Station  5  data  allows  us  to  do  this,  and  we  can  take  into  account,  for 
example,  a  particularly  stormy,  once  in  the  century,  month  that  will 
Introduce  a  profound  change  in  the  delta  similar  to  that  produced  by 
ordering  the  waves.  The  model  should  become  much  more  sophisticated  to 
be  useful  in  practice.  For  example,  the  input  of  river  sediment  as  a 
constant  20,000  m3/day  is  unrealistic,  and  will  vary  with  weather 
conditions  in  the  river  system,  which  may  or  may  not  be  related  to  the 
weather  producing  the  waves  that  is  spreading  the  sand  at  the  delta. 
Further  investigation  is  needed,  which  should  also  include  an  attempt 
to  start  from  an  actual  coastline  instead  of  the  x-axis. 
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APPENDIX  A 

A  COMPUTER  MODEL  OF  DELTA  GROWTH 


List  of  Variables 

Subroutine  INPUT 

KYES  -  constant  filled  with  alphanumeric  value  'YES' 

KCHECK  -  accepts  an  alphanumeric  value  to  compare  with  KYES 

TINTR  -  time  interval  at  which  Delta7  writes  a  coastline,  in  days 

SLOPE  -  the  slope  of  the  beach,  not  in  degrees,  rise  over  run, 

non-dimensional 

CTOTAL  -  the  total  length  of  the  simulation  In  days 

DLIM  -  the  accuracy  of  the  approximation;  when  Ay  is  less  than 

DLIM  the  iteration  stops 

RVRSND  -  the  amount  of  riversand  carried  to  the  delta  in  cubic 

meters/day 

Main  Program 

A  -  is  a  dummy  array  to  write  the  coastline  in 

Y  -  is  the  downdrift  coast  in  meters 

YU  -  is  the  updrift  coast  in  meters 

AB  -  0,  the  angle  between  the  breaking  wave  and  the  coast 

in  radians 

RVRMTH  -  the  angle  of  the  rivermouth  in  radians;  since  the  river 

is  fixed  at  200  meters  width  it's  tangent  is  [ YU( 1 )  -  Y(l)]/200 

DELTAY  -  Ay.  the  accretion  or  erosion  of  a  cell's  height  Y  in  meters 

I  -  an  index 

TEST  -  stores  the  maximum  value  of  AB  in  radians  in  order  to  print 

a  warning  if  it  becomes  unrealistically  large 

T  -  time  in  days 

TOTAL  -  a  counter  that  stores  the  total  time  until  the  next  printed 
coastline  is  due;  in  days 
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J 


an  index 


SANDIN 

SANDOT 

D 

K 


the  inputted  sand  for  a  specific  cell  in  cubic  meters 
per  day 

the  outputted  sand  for  a  specific  cell  in  cubic  meters 
per  day 

the  depth  of  the  cell;  or  the  depth  of  the  sand  which 
is  affected  by  mass  transport, in  meters 

an  index 


TINTR,  SLOPE, 

CTOTAL, 

DLIM  RVRSND  -  same  as  in  subroutine  INPUT 


Subroutine  SWELL 


AO 

P 

TDUR 

L 

KHB 

H 

T 

A 

D 

SINA 

COSA 


-  the  angle  the  refracted  wave  makes  with  the  x-axis 
in  radians 

-  the  power  of  the  wave  in  108  ergs/cm-sec 

-  the  total  length  of  the  simulation  at  which  point  the 
particular  wave  has  passed,  in  days 

-  the  wavelength  of  the  deepwater  wave  in  meters 

-  (kh)b  where  h  is  the  water  depth  and  k  =  2tt/L 

-  deepwater  wave  height,  the  significant  height  in  meters 

-  the  period  of  the  deepwater  wave  in  seconds 

-  a0,  the  angle  the  deepwater  wave  makes  with  the  x-axis, 
input  as  degrees  and  converted  to  radians 

-  the  duration  of  the  wave,  in  days 

-  sin  ab;  used  to  save  an  extra  processing  of  the  sine 
routine 

-  cos  ab;  used  to  save  an  extra  processing  of  the  cosine 
routine 
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APPENDIX  B 

LISTING  OF  SUBROUTINE  INPUT 


C  THIS  IS  A  MODIFIED  DELTA4*  PROGRAM  WHICH  ACCEPTS  LARGE  BLOCKS  OF 
C  WAVE  DATA  TO  CHECK  DCLTA  GROWTH  RATHER  THAN  JUST  ONE  SPECIFIC  UA< 
C  BASE  PROGRAM  DELTA4  IS  STORED  ON  DXOl  AS  DELTA4 .FTN  *  1 
C 
C 

SUBROUTINE  INPUT 

COMMON  T INTR * SLOPE  * CTOTAL * DL IM * RVRSNU 
DATA  K YES/  '  NO ' / 

92  TYPE  1 


1 

FORMAT 

(/*  ' 

AT  WHAT 

INILRVAI.S  DO  YOU  WANT 

PRINTS?' */) 

ACCEPT 
TYPE  2 

101 » 

TINTR 

n 

FORMAT 

(/*  ' 

WHAT  IS 

THE 

beach  slope  ?' »/) 

ACCEPT 
TYPE  3 

201  * 

SLOPE 

3 

FORMAT 

</*  ' 

WHAT  IS 

THE 

SUM-DURATION  OF  THE  SIMULATION?'*./ 

ACCEPT 

101  * 

CTOTAL 

TYPE  22 
ACCEPT 
TYPE  7 

101  * 

RVRSND 

-} 

/ 

FORMAT 

(/*  ' 

WHAT  IS 

THE 

DESIRED  ACCURACY? 

C ABOUT  .011'*/) 

ACCEPT 

701 

*  DL  IM 

TYPE  8 
TYPE  S3 
TYPE  81 

TYPE  82  *  T INTR  *  CTOTAL*  SLOPE*  RVRSND.  DUM 
RVRSND  =  .5  *  RVRSND 
TYPE  9 

ACCEPT  901*  KCKLCK 
IF(KCHECK.EO.KYtS)  GO  TO  92 
TYPE  39 
TYPE  4 
RETURN 

8  FORMAT (/.l OX. 'CHECK  YOUR  INE'U  I  !  !  !  '  ) 

81  F0RMAT<5X, 'DAYS' *15X* 'CUBIC  METERS/DAY') 

02  F0RMAT(2X*2F5.0*4X*F7.3*5X,F8.2*5X*F8.6*/) 

C3  'FORMAT  < '  INTLK  CTOTAL  SLOPE  '* 4X *' RVRSND ', 7X *' DL IM ' ) 
22  FORMAT (/*'  HOW  MUCH  SAND  IN  THE  RIVER  PER  DAY?'*/) 

9  FORMAT ( '  IS  THIS  INPUT  CORRECT?  CYES  OR  NOT'*/) 

901  FORMAT ( A2 ) 

4  FORMAT ( / * 10X * '  DCLTA6  IS  RUNNING*  BC  PATIENT'*/) 

39  FORMAT  < / *  5X  *  '  DX  IS  SET  AT  100  METERS*  DT  A7  ,1  DATS') 

101  FORMAT ( F8 . 2 ) 

201  FORMAT ( F6 . 3 ) 

701  FORMAT  <  F8 . 5 ) 

END 
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APPENDIX  C 

LISTING  OF  MAIN  PROGRAM  DELTA7 


DIMENSION  A ( 53  >  f  Y( 100) »  YU(IOO) 

COMMON  TINTR . SLOPE  rCTOTAL f  DL IM  r R VRSND 
COMMON  /WdVE/AO.P.TDUR 
DOUBLE  PRECISION  AB.RVRMTH.DELTAY 
DO  10.1*1 . 50 
Y ( I  >  *  0. 

10  YU( I )  -  0. 

DO  11  1  =  1 .  53 
J 1  A ( I )  =  0. 

TEST  =  0. 

TDUR  =  0. 

T  =  0. 

CALL  ASSIGN  (2.'WAVC1  OAT',?) 

L 

C  NOTE  THAT  THL  METHOD  0 r  VARYING  WAVE  DATA  IS  TO 
C  CHANGE  THE  CONTENI  OP  THE  WAVE!  DATA  FILE 
C 

CALL  ASSIGN  ( 3 . 'DCLDAT.DAT' >10) 

(.•ALL  SWELL 
CALL  INPUT 
TOTAL  =  TINTR 
C 

100  J  =  1 

RVRMTH  =  <Y( 1 )- YU( 1 ) )*.005 
RVRMTH  =  ATAN ( RVRMTH )  i  AO 
SANDIN  =  RVRSND  -  3425 . *P*S I N ( 2 . *RVRM fH ) 

20  J  =  J  +  1 

AB  =  AT  AN  (  <  Y  <  J  1  )  Y  (  J-i  1  )  )  *  .  005  )  -  AO 

IF  ( A£< .  CT .  TEST  )  TEST® AR 

SANDOT  =  3425.  *  P  *  SIN(2.*AB) 

D  =  4.0  I  SLOPE  *  Y< J  1 > 

DELTAY  =  (SANDIN  -SANDOT >*. 001 /D 
Y(J-l)  =  Y(J-l)  T  DELTAY 
SANDIN  =  SANDOT 

IF (ABS( DELTAY' .GT.DLIM)  GO  TO  20 

C 

200  K  ■  1 

SANDIN  =  RVRSND  •»  3425 .  *P*SIN  (  2  .  *RVRMTH  ) 

220  K  =  K  1  1 

AB  *  ATAN ( ( YU( K  1 )  YU ( K  +  l )  )* . 005  )  +  AO 

IF ( AB . GT . TEST  )  TEST  =  AB 
SANDOT  =  3425.  *  P  *  SIN(2.*AB> 

D  =  4.0  h  SLOPE  *  YU  ( N- 1 ) 

DELTAY  »  (SANDIN  SANDOT )*. 00 1 /D 

YU(K-l)  =  YU(K-l)  +  DELTAY 
SANDIN  =  SANDOT 

IF(ABS(DELTAY)  .GT.DLIM)-  GO  TO  220 
C 

T  *  T  f  .1 

IF ( T . GT . CTOTAL  >G0  TO  20 
If (T.GT.TDUR)  CALL  SWELL 
IF  (T.LT. TOTAL)  GO  10  100 
2S  TYPE  2?  »  T 

29  FORMAT  (/,  5X »  'THE  ' »F6.0.'  DAY  SHORCLINL  PROFILE. 


Y  VALUES 


30 


IF(SLOPE.EQ.O)GU  TO  292 
TYPE  291 i  SLOPE 
GO  TO  293 

291  FORMAT  </*9X»'AT  A  BEACH  SLOPE  OF  ' i.3*/> 

292  TYPE  294 

294  FORMAT(/»9X* 'PLANE  BEACH  ASSUMED*  SLOPE  =  O'*/) 

293  CONTINUE 
TYPE  600 
TYPE  601 

00  30.  I  *  1.  24*  3 

30  TYFE  3ltY(H2)  *  Yil  H>  *  Y(I)  t  YU(I>  *  YU  ( 1*1  )  *  YU(  1+2) 

31  FORMAT  </*  6F8.2) 

IF  <  TEST . GE . 0 . S  >  TYPE  311 

311  FORMAT  < / *  5X  * '  ***4*WARNING*****  AD  MAY  BE  TOO  BIG!!'*/) 
TEST  =  0. 

C 

DO  301*  1=1*25 
A  ( I  128 ) =YU  < I ) 

301  A  < I )  =  Y ( 26  I ) 

WRITE  (3)  <  A ( I ) *  1  =  1  *  53 ) 

Jf  (TOTAL. GE.CTOTADCO  TO  32 
TOTAL  =  TOTAL  +TINTK 
GO  TO  100 

32  STOP 

601  FORMA  T  <  / »  3X  *  '  J  /  2  '  *  6X  *  '  JH  ' » 6X  *  '  J  '  *  7X  *  '  J  '  ►  6X  *  '  J  I  1  '  *  5X  >  '  J-l  2  '  ) 

600  FORMAT ( / » GX *  ' UPDRH  I  '  *  10X*  ' DOWNliRl FT  '  ) 

END 
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APPENDIX  D 

SUBROUTINE  SWELL 


SUBROUTINE  SWELL 
COMMON  /WAVE/AO tF t  TDllff 
REAL*4  L  .  KHB 

READ  <2*111* END=1 10 )  H »  T .  A  .  D 
A  =  .01745  *  A 
F  =  .955  *  T  *  H**2 
TDUR  =  TDUK  »  D 
L  —  1 *56  X 

KHB  =*1.901  *  (CO S(A)  *  <II/L>**2)** 
S I NA  =  KHB  *  SIN(A) 

COSA  =  SORT ( 1 .  -  SINA**2> 

AO  =  ATAN(SINA/C03A) 

110  RETURN 

111  FORMAT ( 4F6 . 1 ) 

END 
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APPENDIX  E 

A  COMPUTER  SIMULATION  OF  DELTA  GROWTH 


Derivation  of  the  approximation  for  ( kh ) ^ 


The  starting  point  Is  that  energy  flux  is  constant  between  wave 
rays  as  the  deepwater  wave  Is  refracted  towards  the  shore  (Komar  (1976)). 


ECnb  =  E  C  n  b  (1) 

0  0  0  0 

where  P  Is  power,  E  Is  energy  as  defined  below,  C  is  wave  group  celerity, 
b  =  distance  between  rays  and 


n 


sinh  kh  J 


(2) 


where  h  is  depth  in  meters  and  k  =  2tt/L,  where  L  is  the  wavelength, 
also  in  meters.  In  deepwater  as  kh  gets  very  large,  nQ  =  1/2,  and 
in  shallow  water,  n  =  1 . 

Komar  (1976)  uses  the  following  expression  for  energy 

E  =  ^  pgH2  (3) 

where  p  is  density,  g  Is  the  acceleration  due  to  gravity  and  H  is 
significant  wave  height.  Using  the  subscript  o  for  deepwater  conditions 
we  may  then  combine  Equations  (1)  and  (3) 

^P9HoVobo  *  pgH2Cnb  (4) 

U 

arranging  terms,  simplifying  and  solving  for  p— 


(5) 
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Refraction  of  the  approaching  waves  over  a  plane  beach  allows  us  to 


write 


cos  a 
o  _  _ o 

b  cos  a 


(6) 


and  thence  for  a  wave  refracting  and  shoaling  on  a  plane  beach: 


r  c  i 

1/2 

'cos  a  " 

0 

0 

L  2Cn  . 

.cos  a  . 

1/2 


(7) 


Now  starting  with  a  general  expression  for  wave  celerity 

C2  *  (g/k)tanh  kh  (8) 

where  all  terms  are  as  defined  above.  We  note  that  in  deepwater,  that 
is  as  h»  1,  tanh  kh  goes  to  one,  e.g.  tanh  6.5  =  0.599996  and 
therefore  for  deepwater 

c0!  -  g/k0  (9) 


Combining  (8)  and  (9) 


[fl  *Ttar,hkh 

(10) 

L  o  J 

Now,  assuming  that  wave  period  stays  constant 

as  the  wavelength 

changes  we  may  write 

r  _  o' 
c  "  IT 

(ID 

L0  k 

(12) 

J°l° 

fl 

Ho7''  o 

(13) 
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where  cr  =  2tt/T  and  k  =  27t/L.  Therefore  we  may  introduce  equations 
(13)  into  (10) 

S-  -  tanh  kh  (14) 

As  we  approach  breaking  conditions,  (kh)b  gets  very  small  (  «  1)  and  we 
may  approximate 


tanh(kh)b  3=  (kh)b 


(15) 


where  the  subscript  b  denotes  breaking.  Combining  Equations  (7),  (14) 
and  (15),  and  remembering  n  =  1  for  shallow  depth 


Hb . 

r 1  i  i 

i  u 

-cos  aQ  " 

Ho  ' 

L2  J 

_  cos  a 

1/2 


06) 


Komar  (1976)  gives  the  following  equation  for  the  wave  height  at 
breaking: 

Hb  =  .142  LQ  tanh2 (kh)b  (17) 

Using  again  the  approximation  of  Equation  (15)  and  dividing  Equation  (17) 
by  Hq  we  have  a  second  relationship  for  Hb/HQ 

-  -142  JjS.  (kh)b2  (18) 

o  no 

Substituting  Equation  (16)  in  (18)  one  obtains: 


1  1 

7  7kh^ 


cos  a0 ' 
cos  ab  _ 


(19) 
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Solving  for  (kh)b  the  following  results: 


(kh) 


5/2 

b 


.142^  L 

o 


1/2 

'  cos  %  1 
.  COS  ab  J 


(20) 


The  term  ab  must  now  be  eliminated  from  Equation  (20)  for  it  to  be 
useful.  We  start  with  Snell's  Law  for  refraction 


s  i  n  a  _ 
C 


sin  a 

o 


(21) 


and  thence  to  the  breaking  condition  which  allows  us  to  use  Equation  (15) 
once  more  to  write 


sin  ab  =  (kh)b  sin  aQ 

which  after  application  of  the  Pythagoreum  Theorem  yields 


(22) 


cos2ab  =  1  -  (kh)h2sin2aQ  (23) 

Raising  Equation  (20)  to  the  fourth  power  we  may  substitute  (23) 


cos2a 

_ o _ 

1  -  (kh),  2  sin2a 
b  o 


(24) 


It  is  difficult  to  solve  for  (kh)b  here,  so  we  introduce  another 
approximation.  Since  sin  aQ  <  1  and  (kh)b  «  1 ,  we  may  neglect  the 
entire  term  (kh)b2  sin2aQ,  and  solve  for  ( kh ) b  on  the  L.H.S.  by  taking 
the  10th  root. 

2/5 

/  H  \  1/5 

(kh)b  =  1 .901  V  /  (cos  ao) 

which  is  the  expression  used  in  this  study. 
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APPENDIX  F 
WAVE  INPUT  DATA 


WAVE 

DATA  HIX 

PRINTED 

PN  19-SEP- 

79  AT 

14:15105  PAGE 

NO 

HEIGHT 

PERIOD 

ANGLE 

DURATION 

- T 

~0/S~~ 

4 , 

~~  0.0 

1.5 

1.5 

2 

0.5 

4. 

2.5 

2.1 

3.6 

- -3 

0 1 5 

"  4v  * 

J  *  V 

2  ♦  u 

6  •  1 

4 

0.5 

4. 

7.5 

2.4 

8.5 

5 

0.5 

4. 

10.0 

1.5 

10.0 

6 

0.5 

5. 

0.0 

1.5 

11.5 

7 

0.5 

5. 

2.5 

2.1 

13.6 

8 

0.5 

5. 

5.0 

2.5 

16 . 1 

9 

0.5 

J  • 

7.5 

2.4 

18.5 

10 

0.5 

5  • 

10.0 

1.5 

20.0 

11 

0.5 

6. 

0.0 

1  .5 

21.5 

12 

0.5 

6. 

2.5 

2.1 

23.6 

13 

0.5 

6. 

5.0 

2.5 

26.1 

14 

0.5 

6. 

7.5 

2.4 

28.5 

15 

0.5 

6. 

10.0 

1 .5 

30.0 

16 

0 . 6 

4  . 

0.0 

1 .5 

31.5 

17 

0.6 

4  . 

2.5 

n  ^ 

23.6 

IS 

0.6 

4. 

5.0 

36. 1 

19 

0 . 6 

4. 

7.5 

2.4 

33.5 

20 

0.6 

4. 

10.0 

1 .5 

40.0 

21 

0.6 

5 « 

0.0 

1 .5 

41 .5 

nn 

0.6 

5  * 

2.5 

2.1 

43.6 

23 

0.6 

5  % 

5.0 

2.5 

46.1 

24 

0.6 

5 » 

7.5 

2.4 

4S.5 

25 

0.6 

5. 

10.0 

1.5 

50.0 

26 

0.6 

6. 

0.0 

1.5 

51 .5 

27 

0.6 

6. 

2.5 

2.1 

52.6 

28 

0.6 

6. 

5.0 

2.5 

56.1 

29 

0.6 

6. 

7.5 

2.4 

5S.5 

30 

0.6 

6. 

10.0 

1.5 

60.0 

31 

0.7 

4. 

0.0 

1.5 

61.5 

37 

0.7 

4. 

2.5 

2,1 

63.6 

33 

0.7 

4. 

5.0 

2.5 

66.1 

34 

0.7 

4  . 

7.5 

2.4 

68.5 

35 

0.7 

4  . 

10.0 

1  .5 

70.0 

36 

0.7 

5. 

0.0 

1.5 

71.5 

37 

0.7  ' 

5. 

2.5 

2.1 

73.6 

38 

0.7 

5. 

5.0 

n  s 

76.1 

39 

0.7 

5. 

7.5 

2.4 

78.5 

40 

0.7 

5  • 

10.0 

1.5 

80.0 

41 

0.7 

6. 

0.0 

1.5 

81.5 

42 

0.7 

6. 

2.5 

2.1 

83.6 

43 

0.7 

A. 

5.0 

2.5 

86.1 

44 

0.7 

6 » 

7.5 

2.4 

88.5 

45 

0.7 

A. 

10.0 

1.5 

90.0 

46 

0.8 

4. 

0.0 

1.5 

91.5 

47 

0.8 

4. 

2.5 

2.1 

93.6 

48 

0.8 

4. 

5.0 

2.5 

96.1 

49 

0.8 

4. 

7.5 

2.4 

98.5 

\JMU4RWi. 


tM.1  O  LI  H  (S  Cl  O  11  I-*  O'  U  O  (.11-*  (Ml  O  LI  O'  LI  O  tn  H  [MU  O  (.1  H  ff-  w  o  Cl  h*  (Kl 
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WAVE 

DATA  MIX 

PRINTED 

ON  1 9 -SEPT¬ 

■79  AT  14J 

15:05  PAGE 

NO 

HEIGHT 

PERIOD 

ANGLE 

DURATION 

50 

0.8 

4. 

10.0 

1.5 

100.0 

51 

0.8 

5. 

0.0 

1.5 

101.5 

52 

0.8 

5. 

2.5 

2.1 

103.6 

53 

0.8 

5. 

5.0 

2.5 

106.1 

..  _54 

..0,3 _ 

.  5, _  7.5, 

108.5 

55 

0.8 

5 . 

10.0 

1.5 

110.0 

56 

0.8  . 

6. 

0.0 

1.5 

111.5 

57 

0.8 

6. 

2.5 

2.1 

113.6 

53 

0.8 

6. 

5.0 

2.5 

116.1 

59 

0.8 

6. 

7.5 

2.4 

118.5 

60 

0.8 

6. 

10.0 

1,5 

120.0 

61 

0.9 

4. 

0.0 

1.5 

121.5 

62 

0.9 

4  . 

2.5 

2.1 

123.6 

63 

0.9 

4. 

5.0 

2.5 

126.1 

64 

0.9 

4  . 

7.5 

2.4 

128.5 

65 

0.9 

4. 

10.0 

1.5 

130.0 

66 

0.9 

5. 

0.0 

1.5 

131.5 

67 

0.9 

5 . 

2.5 

2.1 

133.6 

63 

0.9 

C 

• 

5.0 

2.5 

136.1 

6? 

0.9 

cr 

J  * 

7.5 

2.4 

133.5 

70 

0.9 

5  • 

10.0 

1.5 

140.0 

71 

0.9 

6. 

0.0 

1 .5 

141.5 

72 

0.9 

6. 

n  .  5 

2.1 

143.6 

73 

0.9 

6  . 

5.0 

2.5 

146.1 

74 

0.9 

6 . 

7.5 

2.4 

14E.5 

75 

0.9 

6. 

10.0 

1.5 

150.0 

76 

1.0 

4  . 

0.0 

1.9 

151 .9 

77 

1.0 

4. 

2  •  5 

2.7 

154.6 

78 

1.0 

4. 

5.0 

3.2 

157.8 

79 

1.0 

4  . 

7.5 

3.1 

160.9 

80 

1.0 

4  . 

10.0 

1.9 

162.8 

81 

1.0 

5 . 

0.0 

1.9 

164.7 

82 

1.0 

f 
«J  • 

2.5 

n  “7 
•  / 

167.4 

83 

1.0 

5 « 

5.0 

3.2 

170.6 

84 

1.0 

5 « 

7.5 

3.1 

173.7 

85 

1 .0 

5  # 

10.0 

1.9 

175.6 

86 

1.0 

6. 

0.0 

1.9 

177.5 

87 

1.0 

6. 

2 . 5 

2.7 

180.2 

88 

1.0  . 

6. 

5.0 

3.2 

133. 4 

8? 

1.0 

6. 

7.5 

3.1 

186.5 

90 

1.0 

6. 

10.0 

1.9 

188.4 

91 

1 . 1 

4 . 

0.0 

1.9 

190.3 

92 

1 . 1 

4  . 

2.5 

n  7 

193.0 

93 

1.1 

4. 

5.0 

3.2 

196.2 

94 

1.1 

4. 

7.5 

3.1 

199.3 

95 

1.1 

4. 

10.0 

1.9 

201.2 

96 

1.1 

5. 

0.0 

1.9 

203.1 

97 

1.1 

5  • 

2.5 

2.7 

205.8 

98 

1 . 1 

e; 
w  • 

5.0 

3.2 

209.0 

99 

1 . 1 

5 ♦ 

7.5 

3.1 

212.1 

VJCD  -a  C0O"0  OHH  D.  (.1  O  WH  O'  LIO  U  HOLl  OU  M  CNLIO  C1M 
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WAVE 

DATA  MIX 

PRINTED 

ON  19-SEP 

-79  AT 

14:i5:05  PAGE 

NO 

HEIGHT 

PERIOD 

ANGLE 

DURATION 

100 

1.1 

5. 

10.0 

1.9 

214.0 

101 

1 . 1 

6. 

0.0 

1.9 

215.9 

102 

1.1 

6. 

2.5 

2,7 

218.6 

103 

1.1 

6. 

5.0 

3.2 

221.8 

104 

1 . 1 

6. 

7.5 

3.1 

224.9 

105 

1.1 

6. 

10.0 

1.9 

226.8 

104 

1.2  •  - 

4. 

•  0.0 

1.9 

228.7 

107 

1.2 

4. 

2.5 

2.7 

231.4 

108 

1.2 

4. 

5.0 

3.2 

234.6 

109 

1.2 

4  . 

7.5 

3.1 

237.7 

110 

1.2 

4  . 

10.0 

1.9 

239.6 

111 

1.2 

5. 

0.0 

1.9 

241.5 

112 

1.2 

5  • 

2.5 

2.7 

244.2 

113 

1.2 

5  • 

5.0 

3.2 

247.4 

114 

A  •  A- 

5 . 

7.5 

3.1 

250.5 

115 

1.2 

5 . 

10.0 

1.9 

252.4 

116 

1.2 

6. 

0.0 

1.9 

254.3 

117 

1.2 

UJ  . 

2  •  5 

2.7 

257.0 

113 

1 .2 

6 . 

5.0 

3.2 

260.2 

11? 

1.2 

6. 

-7  rr 

3.1 

263.3 

120 

1 .2 

6. 

10.0 

1.9 

265.2 

121 

1 .3 

4  . 

0.0 

1.9 

267.1 

i 

1.3 

4  . 

n  er 

t  J 

n  “7 

26? .  £ 

123 

1.3 

4. 

5.0 

t  n 

273.0 

124 

1.3 

A  . 

7.5 

3. 1 

276.1 

125 

1.3 

4. 

10.0 

1.9 

278.0 

126 

1.3 

C 

taJ  • 

0.0 

1.9 

"*79  9 

127 

1.3 

r- 
u  « 

2.5 

2.7 

282.6 

128 

1.3 

5 ♦ 

5.0 

3.2 

,OJ.D 

12? 

1.3 

5. 

7.5 

3.1 

">00  o 

130 

1.3 

tr 

♦ 

10.0 

1 .9 

290 . 8 

131 

1.3 

6. 

0.0 

1.9 

292.7 

132 

1.3 

6. 

n  r- 
*_  «  «J 

2 . 7 

295.4 

133 

1.3 

6. 

5.0 

3.2 

298.6 

134 

1.3 

6  . 

7.5 

3.1 

301.7 

135 

1.3 

6. 

10.0 

1.9 

303.6 

136 

1.4 

4  . 

0.0 

1.9 

305.5 

137 

1 . 4 

4. 

2.5 

2.7 

308.2 

138 

1 . 4 

4. 

5.0 

3.2 

311.4 

13? 

1.4 

4. 

7.5 

3.1 

314.5 

140 

1.4 

4.' 

10.0 

1.9 

316 . 4 

141 

1.4 

5 . 

0.0 

1.9 

318.3 

142 

1.4 

5  • 

2.5 

2.7 

321 .0 

143 

1.4 

5 . 

5.0 

3.2 

324.2 

144 

1 . 4 

5  • 

7.5 

3. 1 

327.3 

145 

1.4 

5  • 

10.0 

1.9 

329.2 

146 

i  .  4 

6. 

0.0 

1.9 

331.1 

147 

1.4 

6. 

2.5 

2.7 

333.8 

148 

1.4 

6. 

5.0 

3.2 

337.0 

149 

1.4 

6. 

7.5 

3.1 

340.1 

39 


WAVE 

DATA"  MIX 

PRINTED 

ON  19-SEP- 

79  AT 

14:i5J05  PAGE 

NO 

HEIGHT 

PERIOD 

ANGLE 

DURATION 

150 

1.4 

6  # 

10.0 

1.9 

342.0 

151 

1.5 

4} 

0.0 

0.3 

342.3 

152 

1.5 

4. 

2.5 

0.4 

342.7 

153 

1.5 

4. 

5.0 

0.4 

343.1 

""154 

i  .5 

4  . 

7 .5 

~  6.4 

343.5  ‘ 

155 

1 . 5 

4. 

10.0 

0.3 

343.8 

156 

1.5 

5. 

0.0 

0.3 

344.1 

157 

1.5 

5 . 

2.5 

0.4 

344.5 

158 

1.5 

5. 

5.0 

0.4 

344 . 9 

15? 

1.5 

5  • 

7.5 

0.4 

345.3 

160 

1.5 

5 . 

10.0 

0.3 

345.6 

161 

1.5 

6. 

0.0 

0.3 

345.9 

162 

1.5 

6. 

2.5 

0 . 4 

346 . 2 

163 

1.5 

6. 

5.0 

0 . 4 

346.7 

164 

1.5 

6. 

7.5 

0 . 4 

347. 1 

165 

1 .5 

6  . 

10.0 

0 . 3 

34^.4 

166 

1 .6 

4  . 

0,0 

0.3 

34-'.  7 

167 

1 . 6 

4  . 

O  CT 

0.4 

342.1 

168 

1 .6 

, 

r-  r. 

0 . 4 

348 . 5 

16? 

1.6 

4  . 

7.5 

0.4 

343.9 

170 

1.6 

4  . 

10.0 

0.3 

340.2 

171 

1.6 

5  ♦ 

0.0 

0.3 

340.5 

172 

1 .6 

rr 

2.5 

0.4 

349  .  o 

173 

1  .6 

c 

5.0 

0 . 4 

3Srt .  3 

174 

1 .6 

5 . 

7.5 

0 . 4 

350 . 7 

175 

1 .6 

5 » 

10.0 

0.3 

351 .0 

176 

1 .6 

6. 

0.0 

0.2 

351 . 3 

177 

1.6 

6. 

2.5 

0.4 

351.7 

178 

1.6 

6. 

5.0 

0.4 

352 . 1 

17.9 

1 .6 

6  . 

7.5 

0-4 

7=;  o  e; 

180 

1 .6 

6, 

10.0 

0.3 

352.8 

181 

1.7 

4  . 

0.0 

0.3 

357 , 1 

182 

1.7 

4  . 

2.5 

0.4 

353.5 

183 

1.7 

4. 

5.0 

0.4 

353.9 

184 

1.7 

4  . 

7.5 

0.4 

354.3 

185 

1.7 

4. 

10.0 

0.3 

354.6 

196 

1.7 

er 

*J  ♦ 

0.0 

0.3 

354  .  o 

187 

1.7 

5. 

2.5 

0.4 

355 . 7 

188 

1.7 

5. 

5.0 

0 . 4 

355.7 

189 

1.7 

5 . 

7.5 

0.4 

356.1 

190 

1 . 7 

5 . 

10.0 

0.3 

356.4 

191 

1.7 

6. 

0.0 

0.3 

356, 7 

192 

1.7 

6. 

2.5 

0.4 

357.1 

193 

1.7 

6. 

5.0 

0.4 

357.5 

194 

1.7 

6. 

7.5 

0.4 

357.® 

195 

1  .7 

6 . 

10.0 

0.3 

358.2 

196 

1.8 

4  . 

0.0 

0.3 

358.5 

197 

1.8 

4. 

i  5 

0.4 

75, o  >  o 

190 

1.0 

4  . 

5.0 

0.4 

35?  .  3 

199 

1 .8 

4. 

7.5 

0 . 4 

350.7 

40 


iAVE 

DATA  MIX 

PRINTED 

ON  19-SEP- 

79  AT 

14:15505  PAGE 

NO 

HEIGHT 

PERIOD 

angle" 

DURATION 

200 

1.8 

4. 

10.0 

0.3 

36C.0 

201 

1.8 

5. 

0.0 

0.3 

360.3 

202 

1.8 

5. 

2.5 

0.4 

360.7 

203 

1.8 

5. 

5.0 

0.4 

361.1 

204 

1.8 - 

..  7  . 

n .  a 

361.5 

. . 

205 

1.8 

5. 

10.0 

0.3 

361.8 

206 

1.8 

6. 

0.0 

0.3 

362.1 

207 

1.3 

6. 

2.5 

0.4 

362.5 

205 

1.8 

6. 

5.0 

0.4 

362.9 

20? 

1.8 

6. 

7.5 

0.4 

363  ♦  0 

210 

1.8 

6. 

10.0 

0.3 

363.6 

211 

1 .? 

4. 

0.0 

0.3 

363.9 

n  ^  o 

1.? 

4. 

2  ♦  5 

0,4 

364.3 

213 

1.? 

4. 

5.0 

0.4 

364 . 7 

214 

1  .? 

4  . 

7.5 

0.4 

365.1 

215 

1 .? 

4. 

10.0 

0.3 

365.4 

216 

1.? 

5, 

0.0 

0.3 

6  5  •  / 

217 

1.9 

IT 

U  ♦ 

O  tr 
♦  %J 

0.4 

366 . 1 

212 

1.9 

«J  * 

5 . 0 

0.4 

T  6,  4  tr 

O  O  ♦  w 

219 

1.9 

5  * 

7*5 

0.4 

366 . 9 

'-*20 

1  .? 

«T 

♦ 

10.0 

0.3 

367.2 

221 

1.9 

6. 

0.0 

0 . 3 

-74-7  tr 

900 

1.9 

6. 

2.5 

0.4 

7/0  O 

w *  W-»  /  ♦  7 

223 

1.9 

6. 

5.0 

0.4 

^  w  G  ♦  w 

224 

1.9 

6. 

7.5 

0.4 

368  ♦  / 

°25 

1.9 

6  . 

10.0 

0.3 

369*0 

226 

1.5 

6. 

0.0 

1 . 1 

370 . 1 

02  7 

1.5 

6. 

2.5 

1.7 

371 . 8 

228 

1.5 

6. 

5.0 

1.9 

373.7 

OO^ 

1.5 

6. 

7.5 

1.9 

-7  9«=*  / 

M  .  *J  •  W 

230 

1.5 

6. 

10.0 

1 . 1 

376.7 

231 

1.5 

7. 

0.0 

1 . 1 

377. S 

232 

1.5 

7. 

2.5 

1.7 

379.5 

233 

1.5 

7. 

5.0 

1.9 

381 . 4 

234 

1.5 

7. 

7.5 

1.9 

383.3 

235 

1.5 

7. 

10.0 

1.1 

3S4.4 

236 

1.5 

e. 

0.0 

1 . 1 

385.5 

237 

1.5 

8. 

2.5 

1.7 

387.2 

238 

1.5 

8. 

5.0 

1.9 

38?  .  1 

23? 

1.5 

8. 

7.5 

1.9 

391.0 

240 

1.5 

e. 

10.0 

1  .  1 

392.1 

241 

1.6 

6. 

0.0 

1  .  1 

393 . 2 

242 

1.6 

6. 

2.5 

1.7 

391.9 

243 

1.6 

6. 

5.0 

1.9 

396.8 

244 

1.6 

6  ♦ 

7.5 

1.9 

393.7 

245 

1.6 

6  ♦ 

10.0 

1 . 1 

399.8 

246 

1.6 

7. 

0.0 

1 . 1 

400.9 

247 

1.6 

7. 

2.5 

1.7 

402.6 

248 

1.6 

7. 

5.0 

1.9 

404.5 

24? 

1.6 

7. 

7.5 

1.9 

406 . 4 

41 


WAVE 

DATA  MIX 

PRINTED 

ON  17-SEP- 

79  AT 

14:i5J05  PAGE 

NO 

HEIGHT 

PERIOD 

ANGLE 

DURATION 

250 

1.6 

7. 

10.0 

1.1 

407.5 

251 

1.6 

8. 

1 .1 

408.6 

252 

1.6 

8. 

2.5 

1.7 

410.3 

253 

1.6 

8. 

5.0 

1.9 

412.2 

254 

1.6 

8. 

7.5 

1.9 

414 . 1 

255 

1.6 

8. 

10.0 

1.1 

415.2 

'  256 

~  1.7 - 

-V.”1 - 

•  *  0.0 

“  '1  . 1 

416.3 

257 

1.7 

6. 

2.5 

1.7 

418.0 

258 

1.7 

6. 

5.0 

1.9 

419.9 

259 

1.7 

6. 

7.5 

1.9 

421.8 

260 

1.7 

6. 

10.0 

1  .  1 

422.9 

261 

1.7 

7. 

0.0 

1.1 

424.0 

262 

1.7 

7. 

2.5 

1.7 

425.7 

263 

1.7 

7. 

5.0 

1.9 

427.6 

264 

1 .7 

7. 

7.5 

1.9 

429.5 

265 

1.7 

7. 

10.0 

1 . 1 

430.6 

266 

1.7 

8. 

0.0 

1  . 1 

431.7 

267 

1.7 

8. 

2 . 5 

1.7 

433 . 4 

268 

1.7 

8. 

5.0 

1.9 

435.3 

26° 

1.7 

8. 

7.5 

1.9 

4  37.2 

270 

1.7 

8. 

10.C 

1  . 1 

43G  .  7 

271 

1.8 

6  . 

0.0 

1  .  1 

4  39. 4 

272 

1.8 

L 

’  ♦ 

2.5 

4  “7 

X  ♦  / 

441.1 

-r 

A.  /  U 

1.8 

6. 

5.0 

1.9 

443.0 

274 

1.8 

6. 

7.5 

1.9 

4  4  4  .  c 

275 

1.8 

6. 

10.0 

1  .  1 

446.0 

276 

1.8 

”7 
/  » 

0.0 

1.1 

447.1 

277 

1.8 

•7 
/  * 

O  C" 
m-  ♦  «J 

1 .7 

448 . 8 

278 

1.8 

7. 

5.0 

1.9 

450. ^ 

279 

1.8 

7. 

7.5 

1 .9 

452.6 

280 

1.8 

7. 

10.0 

1  .  1 

453-7 

281 

1.8 

8. 

0.0 

1  .  1 

454.9 

282 

1.8 

8. 

2.5 

1.7 

456 . 5 

283 

1.8 

8. 

5.0 

1 .9 

458 . 4 

284 

1.2 

e. 

7.5 

1 .9 

4  6  o  ,  i 

285 

1.8 

8. 

10.0 

1  . 1 

461.4 

286 

1  .* 

6. 

0.0 

1  .  1 

462.5 

287 

1.7 

6.  . 

2.5 

1.7 

464.2 

288 

1.9 

6. 

5.0 

1.9 

466. 1 

289 

1.9 

6. 

7.5 

1.9 

468 . 0 

290 

1.9 

6. 

10.0 

1  .  1 

4  69  1 

271 

1.9 

7. 

0.0 

1  .  1 

470 . 2 

272 

1.9 

7. 

2.5 

1.7 

471 .9 

273 

1.7 

7. 

5.0 

1.9 

473.2 

294 

1.9 

7. 

7.5 

1.9 

4  75 . 

295 

1.9 

7. 

10.0 

1  .  1 

4  7  6 . 8 

276 

1.9 

8. 

0.0 

1  .  1 

4?~' 

297 

1.9 

8. 

2.5 

1.7 

a  ’’o  ,  £ 

298 

1.9 

8. 

5.0 

1  9 

^91.5 

299 

1  .9 

8. 

7.5 

1.9 

493 . 4 

"*■*  r*C, \  ^  &**-**• 


42 


WAVE 

DATA  MIX 

PRINTED 

ON  19-SEP- 

79  AT 

14115J05  PAGE 

NO 

HEIGHT 

PERIOD 

ANGLE 

DURATION 

300 

1.9 

S. 

10.0 

1 . 1 

484.5 

301 

'  2.0" 

6. 

o.o' 

0.8 

485.3" 

302 

2.0  “ 

'6.  -  ' 

"  2.5 

1.1 

486.4 

303 

2.0 

6. 

5.0 

1.3 

487.7 

304 

2.0 

6. 

7.5 

1.3 

489.0' 

'  303 

2.0 - 

'  6.  - 

10.0  - 

0.8 

489.8 

306 

2.0 

7. 

0.0 

0.8 

490.6 

307 

2.0 

7'. 

2.5 

1  . 1 

491 .7 

308 

2.0 

7. 

5.0 

1.3 

497.0 

309 

2.0 

7. 

7.5 

1.3 

494.3 

310 

2.0 

7. 

10.0 

0.8 

495 . 1 

311 

2.0 

o 

<_/  • 

0.0 

0.9 

495 . 0 

312 

2.0 

8. 

2.5 

1  .  1 

497.0 

313 

2.0 

8. 

5.0 

1.3 

498.3 

314 

2.0 

8, 

7.5 

1.3 

499.6 

315 

2.0 

8. 

10.0 

0.8 

500.4 

316 

2.1 

6. 

0.0 

0.8 

501.2 

317 

2.1 

6  . 

2.5 

1 . 1 

502.3 

318 

2. 1 

6  . 

5 . 0 

1 . 3 

er  -7  / 

319 

2 . 1 

6  . 

7.5 

1 .3 

50^  .  9 

320 

2 . 1 

6. 

10.0 

0.8 

505 . 7 

321 

2.1 

7. 

0.0 

0.8 

-J  t  vJ 

322 

2. 1 

7  . 
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Original  data  had  angles  a0  from  0  to  40°.  These  caused  the 
model  to  "explode"  without  refraction,  i.e.,  for  Delta6. 
Therefore,  all  angles  were  reduced  by  1/4  for  Delta6.  In  order 
to  retain  a  basis  for  comparison  with  the  refracted  wave,  the 
same  angles  were  used  for  Delta7  to  produce  graphs  7,  8,  9  and 
10  (Figs.  9-12).  The  wave  angles  were  returned  to  their  original 
magnitude  for  Graph  11  (Fig.  13). 


INPUT  WAVE  HEIGHTS  -  1ST  QUARTER 


Ordered  wave  data  compared  to  three  random  sorts.  Note  that  there  are 
different  total  durations  for  each  curve  since  each  wave  had  a  different 
duration,  and  waves  were  "shuffled"  with  all  their  parameters  as  one  record 
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Figure  17  Delta7  ordered  increasing  wave  heights.  Wave  height  was  the  primary  method 

of  ordering  waves.  Although  within  height  ranges,  the  waves  were  also  ordered 
by  increasing  period.  Random  wave  data  as  shown  by  the  quarterly  representations 
in  Figures  15  and  16.  Does  not  follow  this  pattern  in  any  way. 


